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Scenarios for controllable creation, trapping and holding of single and multiple solitons in a 
specially designed nonlinear Bragg grating (BG) are proposed. The setting includes a chirped BG 
segment, which is linked via a local defect to a uniform BG with a built-in array of defects. A 
parabolic relation between the trapping position of the incident soliton and its power is obtained. 
Simultaneous trapping of two and three solitons at different locations is demonstrated too. 


OCIS (060. 3735) Fiber Bragg gratings; (060.5530) Pulse 
propagation and temporal solitons; (060.1810) Buffers, 
couplers, routers, switches, and multiplexers 

I. INTRODUCTION 

A possibility to dramatically reduce the veloci^ of ul¬ 
trashort light pulses in various optical materials [l|-Q has 
drawn a great deal of attention, since the observation of 
slow 0 and completely halted light [5[ in experiments. 
However, the controlled storage of this ultrashort light 
pulse has been a challenge, which hinders their use in 
broadband optical signal processing, such as ultrafast 
coding, decoding and multiplexing. Different settings 
were proposed to realized the slow hghb such as elec- 
tromagnetically induced transparency [ll, @, [3 ? simulated 
Brillouin (sf and Raman [9[ scattering, coherent popula¬ 
tion oscillations [H, dispersion engineering in photonic- 
crystal waveguides others. Among these tech¬ 

niques, Bragg-grating (BG) structures, which induce very 
strong dispersion in a vicinity of the photonic bandgap 
M , are used too, for the generation the slow light due to 
its tunability and cascadability in practical applications. 

The relation between the delay and bandwidth of 
pulses was previously experimentally demonstrated in a 
linear resonant systems [I 3 • To support the slow propa¬ 
gation of non-spreading pulses in the BG, the Kerr non¬ 
linearity may be used, as it is not subject to the limitation 
of the trade-off between the pulse delay and bandwidth, 
which is inherent in the linear settings [12] . Standing and 
moving robust light pulses (BG solitons), whose velocity 
may be orders of magnitude smaller than the velocity of 
light in the host material, were studied in detail theoreti¬ 
cally [iM3- Producing quiescent pulses by collisions of 
counterpropagating BG solitons was predicted too. 

In experiments, BG solitons, with the velocity corre¬ 
sponding to 76% of the speed of li ght in the uniform 
medium, were first produced in Ref. [^. Subsequently, 
a number of experiments [niliiil were carried out to 
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investigate slow-light effects in BG structures. In par¬ 
ticular, high intensity of the pulse leads to a change 
of the local refractive index via the Kerr nonlinearity, 
which helps to enhance the retardation of pulses prop¬ 
agating in the BG structure. Recently, large femtosec¬ 
ond Kerr nonlinearity was ex per imentally demonstrated 
in cholesteric liquid crystals Its large Kerr coeffi¬ 

cient would considerably reduce the high intensity. The 
lowest experimentally demonstrated velocity of solitons 
in BGs is O.I6co [ll|, where cq is the speed of light in 
vacuum. 

However, there remain challenging problems impeding 
the use of the slow light pulses in BGs. First, the creation 
of very slow or standing optical modes in BG structures 
has not yet been reported in experiments. Second, to the 
best of our knowledge, no work has demonstrated, as yet, 
controllable trapping and storage of single and multiple 
pulses at different positions in BGs, a reason being that 
it is difficult to create suitable input pulses for this pur¬ 
pose, using uniform BGs. Indeed, as the optical fields 
inside the BG include forward- and backward-traveling 
waves, generating such very slow or quiescent pulses re¬ 
quires to couple forward and backward waves with nearly 
equal powers into the BG waveguide, doing which for the 
backward wave being obviously difficult. 

Recently, we have proposed a specially tailored struc¬ 
ture, built of a linearly chirped BG segment linked, 
through a local defect, to a uniform grating [26|. This 
setting demonstrates a possibility to create very slow or 
standing stable pulses, the advantage being that the pulse 
coupled into the uniform grating may be manipulated by 
means of the chirped segment. This scheme makes it pos¬ 
sible to prepare the right mix of forward- and backward¬ 
traveling fields. Using it, the creation of BG solitons 
with extremely small and zero velocities was predicted, 
offering an essential improvement of previously published 
results [ 23 , [ 13 . 

In this work, we aim to develop the above-mentioned 
setting by introducing a periodic set of defects into the 
uniform grating, to achieve controllable trapping and 
storage of stable single and multiple pulses at different 
positions in the BG. In section H, we present the nec¬ 
essary model, which is based on nonlinear coupled-mode 
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equations (CMEs) for the BG. In section III, the oper¬ 
ation of the system is demonstrated by means of CME 
simulations. The paper is concluded by Section IV. 

II. THE MODEL 

The purpose of introducing a chirped BG before the 
uniform grating with a periodic array of defects is to pro¬ 
vide a suitable setting for the generation of slow pulses, 
which can be subsequently trapped in the main section of 
the system. Eirst, the setting without the defect lattice 
is defined by means of the spatially modulated profile of 
the local refractive index [^: no[l + 2An cos(27r2;(l + 
Cz)/K q)] for (0 < z < Lc), and no[l + 2Ancos(27rz(l -h 
d/c)/Ao)] for (Lc < z < L). Here z is the propagation 
distance, L the total length of the setting, Lc the length 
of the chirped Bragg gratings, and C the chirp coefficient. 
Eurther, Aq is the BG period at the input edge, no the 
average refractive index, and An the amplitude of the 
refractive-index modulation. 

The nonlinear CME system is an accurate model for 
the propagation of optical pulses in BGs with the Kerr 
nonlinearity and inhomogeneity of the refractive index 

[Hiiii: 



FIG. 1. (Color online) A schematic of the system, built of the 
linearly chirped BG segment (on the left-hand side) followed 
by the uniform grating with an inserted periodic array of local 
defects, which is described by the set of parameters (S', dw^s). 


balance between the powers of the two waves, which is 
necessary for the creation of a very slow BG soliton in the 
uniform-BG [iM3- The soliton’s velocity is determined 
by the remaining imbalance between the powers: 


±i 


dEf,t 

dz 


i dEf^i) 

Vg dt 


-h 6{z)Ef^i) -h k,{z)Ei)J 


+7(l-E’/.b|^ + 2|-E’b./|^)-E’/,b = 0, 


( 1 ) 


where t is time, Ef and are amplitudes of the 
forward-traveling and backward-traveling waves, respec¬ 
tively, Vg = cq! no is the group velocity in the material of 
which the BG is fabricated, and 7 = n 2 Co’/c is the nonlin¬ 
earity coefficient, with uj being the carrier frequency and 
77-2 the Kerr coefficient. The strength of the coupling be¬ 
tween forward and backward waves is hc = 7rAn/Ao, while 
S{z), which represents the local wavenumber detuning, is 

. _ f ^0 “ 27rC2;/Ao, for 0 < z < 

\ (5o — 27rCI/c/Ao, for Lc < z < L, (2) 

where (5o = 27rno/A — tt/Aq is the detuning at the left 
end of chirped segment. The nonlinear terms in Eqs. 
m account for the self- and cross-phase modulation, re¬ 
spectively. We denote the chirped segment of length Lc 
as chirped-BG (the chirped Bragg grating), and the uni¬ 
form one as uniform-BG (the uniform Bragg grating), 
the length of which is L — Lc. The setting implies that 
the local BG period decreases with z in the chirped-BG, 
hence the nominal Bragg wavelength is gradually shifted 
to smaller values for pulses running through the chirped 
grating. In other words, the wavelength of the input 
pulse, originally taken near the blue edge of the bandgap, 
drifts into the depth of the bandgap in the course of the 
propagation of the pulse [26|. As a consequence, the 
power of the forward-propagating wave is gradually con¬ 
verted into the backward wave, helping to introduce the 


Vsol=Vg{l-f)/{l + f) (3) 


fo lEf(z,t)l2dz 

By introducing a defect at the junction between the 
chirped-BG and uniform-BG segments, one can create a 
standing soliton stably trapped by the defect [26[. Here, 
we introduce an array of defects in the uniform-BG, as 
multiple positions at which the stable trapping should be 
possible. Eor this purpose, we define the local detuning 
in the uniform-BG region as 


Sd{z) = 5{z)x 

N . 

1 + e y] exp p 

n=l ^ 

[z-{l + {n-l)S)]^\ 

dl I 



(4) 


with 5 being the depth of each defect, N the total num¬ 
ber of defects, S the spacing between them, and the 
width of defect. The array is displayed in Eig. [U which 
shows 5(i{z) as a function of z. Thus, the defect array is 
determined by a parameter set, (E, dw^e). The respective 
shape of the effective pot ential for solitons, V(z), which is 
opposite to Sd(z) [ 31 I, l32|, predicts that, when the pulse’s 
kinetic energy, E, is smaller than the local height of the 
effective potential, Vb, the pulse cannot pass the local po¬ 
tential maximum. If K slightly exceeds Vb, it is expected 
that the pulse may pass the maximum into the adjacent 
potential well at a small velocity, and get trapped there 
1^ . Below, this possibility is confirmed by systematic 
simulations of Eqs. m- 
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III. NUMERICAL ANALYSIS 
A. Controllable trapping of the single pulse 

In this section, systematic simulations of the model are 
reported for physical parameters of silicon, whose Kerr 
coefficient and average refractive index are 4.5 x 10“^^ 
cm^/W and no = 3.42 [33|. The amplitude of the 
refractive index-modulation is taken as An = 0.006, 
a typical chirp coefficient in the chirped-BG segments 
C = 2.5064 X 10“^ cm“^, and the BG period at the in¬ 
put edge of the sample is fixed to be Aq = 154.1 nm. 
Simulations based on Eqs. (HD and dH) were carried out 
with the boundary condition corresponding to the Gaus¬ 
sian pulse. 


Ef{z = 0,t) = Eo exp [-{t/to)'^] , (5) 



with temporal width to = 16 ps, launched into the sys¬ 
tem at carrier wavelength 1053 nm, which is placed near 
the blue edge of the photonic bandgap. Further, the total 
length of the system and the length of the chirped-BG 
segment are fixed as L = 1.54 cm and Lc = 0.193 cm, 
respectively (below, the influence of Lc on the pulse stor¬ 
age is discussed), thus the length of the uniform-BG seg¬ 
ment is 1.347 cm. The number of defects in the array is 
A = 10, with individual defect characterized by param¬ 
eters {S^dw^e) = (0.132 cm, 50 //m, 0.04), i.e., the space 
between the defects, S', is much larger than the width of 
each one, d^. In comparison with the spatial width of 
the pulse (the temporal duration of 16 ps corresponds to 
4.8 mm in space) the defect of width 50 jam may indeed 
be considered as a point-like object. 

Simulations of Eqs. (HD were performed, by means of 
the fourth-order Runge-Kutta method. The intensities 
Ip used in simulations are below the damage thresh¬ 
old [33 [Ip is the peak intensity of pulse, defined as 
Ip = \Eo\^, see Eq. ( 0 )]. First, setting the initial pulse 
intensity Ip to 2.07 GW/cm^, trapping of a stable pulse 
by the first defect (i.e., the formation of a standing opti¬ 
cal soliton) is observed, as shown in Figs. [2ja,b). These 
results demonstrate that the Kerr nonlinearity in the 
BG, at the present (quite realistic) intensity level not 
only compensates the pulse’s dispersion in the vicinity 
of the edge of photonic bandgap, securing the formation 
of the gap soliton, but also shifts the photonic bandgap 
to a longer wavelength, which enables the pulse to pass 
the chirped-BG segment, before getting trapped in the 
uniform-BG section. Indeed, when the pulse reaches 
the junction between chirped-BG and uniform-BG, the 
power of the forward and backward components approx¬ 
imately equal, giving rise to a nearly halted pulse, as seen 
from Eq. (ED, which may be readily captured by the first 
potential well. 

The evolution of the pulse’s overall velocity, Vsop de¬ 
fined according to Eq. (ED, is additionally shown in the 
inset of Fig. EJa). It is seen that Vsoi rapidly drops at 
t < 3 ns, due to the transfer from forward component to 




FIG. 2. (Color online) The nonlinear propagation of pulses 
in the system with the defect array with {S,dw,£) = 

(0.132 cm, 50 /xm, 0.04). The first and second columns demon¬ 
strate, severally, the intensity evolutions and the correspond¬ 
ing profiles at t = 12 ns, for three different injection inten¬ 
sities: (a,b) 2.07 GW/cm^ (c,d) 2.72 GW/cm^ (e,f) 3.36 
GW/cm^. Insets of (a), (c), (e) represent the evolutions of 
the pulse’s velocity Vsoh calculated by Eq. ©. 


the backward one in the course of passing the chirped-BG 
segment. After that, Vsoi approaches zero, suggesting the 
trapped of the soliton in the potential well. 

With the incident-pulse intensity increased to Ip = 
2.72 GW/cm^, trapping of a stable soliton in the second 
potential well is observed in Figs. [2fc,d). These results 
reveal that the pulse with this level of the intensity passes 
the first defect, and then keep moving at a small velocity, 
eventually being trapped by the second well. The further 
increase of the intensity to Ip =3.36 GW/cm^, as shown 
in Figs. leads to the trapping of a stable soliton 

in the sixth potential well. 

The above results suggest that the number of the trap¬ 
ping effective potential well increases with the increase 
of the incident-pulse’s intensity. To clearly demonstrate 
this trend, simulations were carried out with different in¬ 
tensities Ip, varying from 2 GW/cm^ to 4.05 GW/cm^, 
for two different defect strengths, 5 = 0.04 and 5 = 0.06. 
The outcome is depicted in Fig. [3l^a), which shows the 
trapping position, n, as a function of the intensity. In 
both cases, the relation between n and Ip can be ap- 
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FIG. 3. (Color online) (a) The relation between the trapping 
position and the initial pulse’s intensity Ip for (S,dw,s) = 
(0.132 cm, 50 /xm, 0.04) and (0.132 cm, 50 /xm, 0.06) (blue and 
red curves, respectively), (b) The ratio of the residual total 
power of the trapped pulse to the input power at e = 0.04, 
versus the initial intensity. 


proximated by 


n = l+a{Ip - Icrf, (6) 


where Icr is a critical value of the injected intensity, Ip^ 
for trapping the pulse in the first potential well, and 
a is obtained by fitting the numerical data. For in¬ 
stance, a = 2.55 cm^/(GW)^ for 5 = 0.04, and a = 1.23 
cm^/(GW)^ for 5 = 0.06. Although the critical in¬ 
tensities are practically identical in both these cases, 
Icr = 2.07 GW/cm^, see Fig. [3Ka), there is a trend to 
decrease of Icr with the decrease of 5. Indeed, in the limit 
of 5 ^ 0 any pulse is able to pass the defect. Naturally, 
pulses stay trapped in a particular potential well when 
Ip takes values in some interval [^. For example, at 
5 = 0.06, the pulse is trapped by the first well for 2.07 
GW/cm^ < Ip <2.72 GW/cm‘^. 

The quadratic dependence displayed by Eq. m can be 
explained, making use of the soliton’s equation of motion. 
In the framework of the CME system (IT]), the soliton’s 
momentum is defined as (see, e.g.. Ref. |1^ ) 


M = i 


‘Ci 


'9E} 

dz 


■Ef 




dz. 


(7) 


As follows from the boundary conditions (|5j), the initial 
momentum of the injected pulse scales as Mq ^ Ip. On 
the other hand, for a relatively slow soliton, its velocity 
is proportional to the momentum, Vsoi ^ Af, see Eq. m 
[l^ . Eurther, the radiative braking force acting on the 
soliton, Fbr, can be estimated as a loss of the momentum 
by the soliton passing a single defect, AM, times the 
number of defects passed in a unit of time, ^ Vsoi/S^ 
see Eq. 0. The perturbation theory [35| shows that 
AM is proportional to the time necessary for the soliton 
to pass the defect, i.e., AM ^ dwjVsoh Thus, we obtain 
Fhr ^ dw /which does not depend on velocity VsoU 
the corresponding equation of motion, dM/dt = —Fbr, 
predicts that the soliton will come to a halt, M = 0, 
at time thait r\j Mq. Einally, an elementary mechanical 
analysis demonstrates that the total distance passed by 
the soliton under the action of the constant braking force 


is ^hait ^ Mq /Fbr ^ Ip , which explains the quadratic fit 
provided by Eq. 

The soliton relaxing to the eventually trapped station¬ 
ary state loses a part of its power through emission of 
linear waves (radiation). Accordingly, the final value of 
the total power of the trapped optical pulse. 



\Ef[z^t oc)| + \E}y(z^t oo)| 


dz 


( 8 ) 


(the integral is computed over a region where the soli¬ 
ton is eventually trapped), normalized to the total input 
power, P(0), is displayed, as a function of the input in¬ 
tensity, /p, in Eig. [3]^b). A trend of Pt/P(0) to decrease 
with the increase of Ip is explained by the fact that the 
pulse with larger power passes a larger number of defects, 
hence it is subject to stronger radiation losses. This curve 
also shows a weak oscillatory behavior with the increase 
of Ip, ratio Pt/P(0) staying nearly constant in narrow in¬ 
tervals of Ip. This phenomenon can be explained: pulses 
with Ip taking values in some intervals, which indeed 
are narrow, stay eventually trapped in the same poten¬ 
tial well. Past these intervals, the ratio Pt/P(0) again 
decreases with the increase of Ip. 

The formation of the gap soliton from the initial pulse 
is additionally illustrated in Eig. |4]^a) by a typical ex¬ 
ample of the temporal evolution of its effective squared 
spatial width, defined as 


lo {\Ef{z,tp + \Eb dz 

/cf {\Ef{z,tp +\Eb{z,tp] dz 


The width relaxes to a constant value by t = 4 ns, when 
the pulse is already trapped by a local potential well, 
and the soliton is well formed. Here, an estimate for the 
distance travelled by the pulse before it forms a trapped 
soliton is presented. Because Vsoi oscillates around zero 
before the formation of a soliton completes [see the in¬ 
set in Eig. He)], we assume that the average velocity of 
the pulse is 0.015c, which corresponds to an estimate for 
the travelled distance ^1.8 cm, which is far larger than 
the effective nonlinearity length, Lnl = 1/ (t^p) — 0.14 
cm. Length Lc of the chirped-BG segment strongly af¬ 
fects the subsequent pulse trapping in the uniform-BG 
section. To illustrate this feature, a relation between 
the trapping site, n, and Lc is displayed in Eig. |4][b), 
for fixed parameters (S', d^^e) = (0.132 cm, 50 /xm, 0.04), 
and fixed incident intensity. Ip = 2.72 GW/cm^. In this 
case, the trapping occurs when Lc takes values between 
0.16 cm and 0.21 cm, with the trapping position, n, de¬ 
creasing with the increase of Lc in this interval. This is 
explained by the fact that the increase of Lc leads to a 
reduction to the pulse’s velocity, hence it is captured ear¬ 
lier, at smaller n. However, the pulse cannot be trapped 
at Lc < 0.16 cm, and at Lc > 0.21 cm. In the for¬ 
mer case, the pulse’s velocity remains too large, allowing 
it to pass the entire uniform-BG segment without being 
halted; in the latter case, the pulse cannot reach the junc¬ 
tion between the chirped-BG and uniform-BG segments. 
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(a) (b) (a) (b) 


FIG. 4. (Color online) (a) The temporal evolution of the 
effective squared width of the pulse, defined as Eq. (|9|), with 
the same parameters as in Figs. [2Uc,d). (b) The trapping 
site, n, as a function of length Lc of the chirped-BG segment, 
for the incident intensity /p = 2.72 GW/cm^ in the structure 
with {S,dw,£) = (0.132 cm, 50 /xm, 0.04). 


eventually suffering decay in the chirped-BG segment, 
due to the strong BG-induced dispersion. 


B. Multi-pulse trapping 

The present BG system is capable to trap several pulses 
too, in different potential wells, provided that the pulses 
are launched into the system successively. For instance, 
two-pulse trapping was demonstrated for the same sets of 
parameters as considered above for the single-pulse case, 
i.e., {S^dw^e) = (0.132 cm, 50 /xm, 0.04). Specifically, 
the first pulse with intensity Ip = 2.72 GW/cm^ [it is 
expected to be trapped at the second potential well, ac¬ 
cording to Figs. [2](c,d)], is launched into the system at 
t = 0, which is followed by launching the second pulse 
with Ip = 2.07 GW/cm^ at t = 3 ns [Figs. [2](a,b) sug¬ 
gests that it may be captured by the first potential well]. 
The outcome of the simulations demonstrates, in Figs. 
5(a,b), that these two pulses are indeed stably trapped 
in two different wells. 

An example of three-pulse trapping is shown in Figs. 
[5](c, d). The first pulse, with Ip = 3.04 GW/cm^, enter¬ 
ing the chirped-BG input edge at t = 0, is trapped by the 
third wells. The second pulse with Ip = 2.78 GW/cm^ 
is launched at t = 3 ns, and is captured by the second 
well. The third pulse with Ip = 2.07 GW/cm^, launched 
at t = 6 ns, is eventually trapped in the first well, as 
suggested by the above results for single pulses. 

Although the above examples of the two- and three- 
pulse trapping seem as simple superpositions of the 
single-pulse dynamical processes considered above, it 
should be stressed that, loosing a considerable part of 
its power on its way to the final halt [see Fig. [3][b)], 
each pulse leaves a trace of emitted radiation in its wake. 
Nevertheless, Fig. [5] demonstrates that the passage of the 
areas “contaminated” by the radiation does not perturb 
the propagation of the secondary solitons, which addi¬ 
tionally attests to the solitons’ stability. 

We also examined the case when the second pulse was 
launched right after the first one, before the first pulse 



FIG. 5. (Golor online) Multi-pulse trapping in the Bragg- 
grating structure with {S,dw,s) = (0.132 cm, 50 /xm, 0.04). 
(a) The trapping of two pulses with incident intensities of 
the first and second pulses Ip — 2.72 GW/cm^ and 2.07 
GW/cm^. The second pulse is launched into the gratings 
at t = 3 ns. (c) The trapping of three pulses, with intensities 
Ip = 3.04 GW/cm^, 2.78 GW/cm^, and 2.07 GW/cm^. The 
pulses are launched into the gratings at t = 0, t = 3 ns, and 
t = 6 ns, respectively. (b,d) Local-power profiles for both sets 
at t = 12 ns. 



FIG. 6. (Golor online) The two-pulse interaction in the system 
with {S,dw,s) = (0.132 cm, 50 /xm, 0.04). (a-c) The incident 
intensities of the first and second pulses /p = 2.72 GW/cm^ 
and 2.07 GW/cm^, while in (d-f) they are /p = 2.07 GW/cm^ 
and 2.72 GW/cm^, respectively. The second pulse is injected 
at t = 0.25 ns (a,d), 0.75 ns (b,e), and 1.0 ns (c,f). 


gets trapped, as in the example shown in Figs. EJa-c). 
In this case, it is concluded that, when the time inter¬ 
val between these incident pulses is < 1 ns, the second 
pulse bounces back, see Figs. [6][a,b). On the other hand, 
if the time interval between the pulses exceeds 1 ns (al¬ 
lowing the first pulse to get well trapped), see Fig. [6][c) 
[also Fig. [Sfa)], these pulses can end up being stably 
trapped in different potential wells simultaneously. In 
addition, the situation where the incident intensities of 
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FIG. 7. (Color online) The two-pulse interaction in the system 
with the reduced spacing of the lattice defect: (S,dw,s) = 
(0.044 cm, 50 /xm, 0.04) (a); (0.066 cm, 50 /xm, 0.04) (b); 
(0.088 cm, 50 /xm, 0.04) (c). The first and second pulses are 
launched at t = 0 and t = 4.1 ns, with initial intensities 
Ip = 2.72 GW/cm^, and /p = 2.07 GW/cm^, respectively. 

the first and second pulses are Ip = 2.07 GW/cm^ and 
2.72 GW/cm^ is demonstrated in Figs. [6l^d-f). In this 
case, the two pulses collide coherently in the course of the 
evolution. Figure [6l^d) demonstrates that these pulses 
interact attractively, finally being trapped in the same 
potential well. On the contrary. Figs. [6l^a,b,e,f), reveal 
repulsive interaction. Moreover, in Fig. Ele), the repul¬ 
sive interaction actually leads to release of an initially 
trapped soliton, which may be used for retrieval of data 
bits stored in the system, in terms of applications. The 
attractive and repulsive sign of the interaction is well ex¬ 
plained by the relative phase, namely, they attract or 
repel each other when they are in phase or out of phase, 
respectively (sH, [36| . 

It is possible to introduce interaction between solitons 
trapped by adjacent defects, reducing the spacing. S', of 
the defect lattice [see Eq. ©]• Denser lattices are of 
obvious interest too in term of applications. In the Figs. 
[71[a,b), and (c), reduced S is taken as S = 0.044 cm, 
S = 0.066 cm and, additionally, S = 0.088 cm, in com¬ 
parison with S = 0.132 cm in Figs. [5] and [6l while 
and 5 keep the same values as before. Although this 
reduction of the spacing is not dramatic, it is sufficient 
to allow actual interactions between solitons originally 
trapped in close potential wells, due to significant over¬ 
lapping between solitons’ tails. For these settings, sim¬ 
ulations were again performed for the first pulse with 
Ip = 2.72 GW/cm^, injected into the system at t = 0, 
and the second pulse with Ip = 2.07 GW/cm^, injected 
at t = 4.1 ns. The large temporal delay allowed the first 
pulse to settle down into a trapped state before the sec¬ 
ond one would appear in the vicinity. Figures [71[a) and 
(b) reveal that the interaction between solitons, which 


are trapped, for a short time, by close defects, indeed 
occurs in these cases, and changes the results. 

In the case of S = 0.044 cm, shown in Fig. [7|[a), the 
second pulse is originally trapped at distance 2S from the 
first one. In this case, the attractively interacting pulses 
hop together into a higher potential well and eventually 
merge in it. For S = 0.066 cm. Fig. [7|[b) demonstrates 
that the first soliton, originally trapped in the second 
well, is pulled by the incident second pulse back to the 
first well, where they merge into a single pulse staying in 
the first well. On the other hand. Fig. [7l[c), correspond¬ 
ing to S = 0.088 cm, shows that the first pulse stays 
trapped in its original position, despite the repulsive in¬ 
teraction with the second incident pulse. As mentioned 
above, the attractive [Figs. [7l[a,b)] or repulsive [Fig. [7l[c)] 
sign of the interaction is determined by the relative phase 
of the two pulses [H, . 


IV. CONCLUSIONS 

We have introduced a system engineered as a concate¬ 
nation of linearly chirped and uniform BGs (Bragg grat¬ 
ings), with the array of local defects embedded into the 
uniform grating. The system of CMEs (coupled-mode 
equations) have been used to simulate the evolution of 
single and multiple pulses injected into the system. The 
systematic analysis has demonstrated that, selecting pa¬ 
rameters of the systems and the intensities of incident 
pulse, the conversion of the pulse into a well-formed BG 
soliton and, eventually, its trapping at a desired position, 
by one of the local defects, are provided by the system. 
The relation of the trapping position and input inten¬ 
sity was found, being close to a quadratic form. The 
co-trapping of two or several pulses may be complicated 
by the strong interaction between them, in the case when 
the temporal delay between the pulses, or the spacing be¬ 
tween adjacent local defects, is relatively small. 

As a further development of the analysis, it may be in¬ 
teresting to study in detail release of a trapped soliton by 
an incident one, an example of which is displayed in Eig. 
He). It may also be interesting to consider an effects of 
a frequency shift added to the input pulse, which corre¬ 
sponds to multiplying input (6) by exp{—iujt) ^ with con¬ 
stant frequency uj. On the other hand, for very slowsoli- 
tons, taking into account optoacoustic effects mediated 
by electrostriction labial may improve the accuracy of 
the model. 
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